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Abstract. We apply the Direct Simulation Monte Carlo (DSMC) method, 
developed originally to calculate rarefied gas dynamical problems, to 
study the gas flow in an accretion disc in a close binary system. The 
method involves viscosity and thermal conduction automatically, and can 
thus simulate turbulent viscosity (a viscosity). 



1. Introduction 

The Direct Simulation Monte Carlo method was originally developed to treat rar- 
efied gas flow (see Bird 1994; Matsuda, Mizutani & Boffin 2001). In DSMC, mi- 
croscopic molecular motions and their mutual collisions are calculated. Macro- 
scopic quantities such as density, velocity and temperature of fluids are cal- 
culated by taking some mean of microscopic quantities. In DSMC, molecular 
collisions are handled in a stochastic manner, hence the Monte Carlo denomina- 
tion of the method. The collision model in DSMC is such that the molecules are 
assumed to be hard billiard balls. Since we solve the molecular motions, molecu- 
lar viscosity and molecular thermal conduction are automatically included in the 
system. This 'physical' viscosity may mimic a turbulent viscosity (a viscosity) 
assumed in the standard theory of accretion disk. 



2. Application of DSMC to a three-dimensional accretion disc 

We consider a binary system with mass ratio q = 0.5. We work in the rotational 
frame. We restrict our computational region to a square, the center of which 
is the primary and just including the LI point from which gas particles are 
injected constantly into the computational region. The region is divided into 
100 x 100 x 20 cells. The cell size is such that dx = dy = 2 dz. We start 
injection of gas particles at a constant rate per time step and follow the time 
evolution of the system. The number of particles in the computational region 
increases with time and reaches a steady state. After a steady state is reached, 
we further calculate to take a time average. Typically we computed up to t=60- 
100 in dimensionless unit, that is about 10-15 orbital periods. 
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Figure 1. Iso-density surface of an accretion disc with velocity vec- 
tors on the rotational plane; (left) 6 particles are injected every time 
step; (right) 30 particles are injected every time step 



In accretion discs, radiative cooling must occur, and some mechanism of 
cooling is essential in our calculation. We here consider a simple model of 
cooling: we reduce the thermal speed of the molecules by a factor / = 0.97 
at each time step. 

3. Results 

Figure 1 displays the time averaged iso-density surface, on which temperature 
distribution is pasted, and velocity vectors in the rotational plane of our sim- 
ulation. Two cases of injected particles number, rii n j = 6 and rii n j = 30, are 
compared. We can observe the gas stream from the LI point, the accretion disc 
and spiral-shaped shocks. In the case of rii n j = 6, we find that the accretion disc 
is very circular. The LI stream hitting at the disc is deflected from its direction 
downward by the disc flow. Because of the collision between the LI stream and 
the disc flow, the iso-density surface is elevated and the temperature on the ele- 
vated area is raised (a hot line or heated wall). In the case of rii n j = 30, the LI 
stream is so strong that the direction of the disc flow is deflected and the disc 
shape is elongated. We may observe a pair of spiral shocks in this case. 
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